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Abstract 

We consider a diffuse interface model for tumor growth recently proposed in [3]. In this new 
approach sharp interfaces are replaced by narrow transition layers arising due to adhesive forces 
among the cell species. Hence, a continuum thermodynamically consistent model is introduced. 

The resulting PDE system couples four different types of equations: a Cahn-Hilliard type equation 
for the tumor cells (which include proliferating and dead cells), a Darcy law for the tissue velocity 
field, whose divergence may be different from 0 and depend on the other variables, a transport 
equation for the proliferating (viable) tumor cells, and a quasi-static reaction diffusion equation 
for the nutrient concentration. We establish existence of weak solutions for the PDE system 
coupled with suitable initial and boundary conditions. In particular, the proliferation function at 
the boundary is supposed to be nonnegative on the set where the velocity u satisfies u ■ > 0, 

where v is the outer normal to the boundary of the domain. We also study a singular limit as the 
diffuse interface coefficient tends to zero. 

Key words: tumor growth, diffuse interface model, Cahn-Hilliard equation, reaction-diffusion 
equation, Darcy law, existence of weak solutions, singular limits. 
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1 Introduction 

Mathematical modeling and analysis of tumor growth processes give important insights on cancer 
growth progression. The models are expected to help to provide optimal treatment strategies. The 
behavior of tumors is a complex biological phenomenon, influenced by many factors, such as cell- 
cell and cell-matrix adhesion, mechanical stress, cell motility and transport of oxygen, nutrients and 
growth factors. In recent years, many mathematical models of cancer have been proposed and various 
numerical simulations have been carried out (cf., e.g., the recent reviews [U [3 |H1 E] ). A variety of 
models are available to investigate different characteristics of cancer: single-phase continuum and 
multiphase mixture models, and methods that combine both continuum and discrete components (cf., 
e.g., [3 Chap. 7]). 

We will address the problem of existence of weak solutions for a PDE system for a tumor 
growth model introduced in [3] (cf. also m and [HD and analyze a singular limit of that model. The 
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works listed above can be framed in the continuum tumor growth models category. This modeling 
approach has become central in the studies of tumor development in applied mathematics (cf. also 
Din])- Actually, the translation of biological processes into models generally turns out to be simpler 
for discrete models than for continuum approaches. Nevertheless, discrete models can be difficult to 
study analytically because the associated computational cost rapidly increases with the number of 
cells modeled. This makes it difficult to simulate millimeter or greater sized tumors. For this reason, 
in larger scale systems (millimeter to centimeter scale), continuum methods provide a good modeling 
alternative. Mixture models, on the other hand, provide the capability of simulating in detail the 
interactions among multiple cell species. 

In the framework of continuum models, the diffuse interface method turns out to be particu¬ 
larly useful to describe multi-species tumor growth processes. In this approach the sharp interfaces 
are replaced by narrow transition layers arising due to the adhesion forces among different cell-species. 
This choice is quite effective since it avoids to introduce complicated boundary conditions across the 
tumor/host tissue and other species/species interfaces. This would have been the case when consid¬ 
ering sharp interface models. Moreover, the diffuse interface approach eliminates the need of tracking 
the position of the interfaces, which is one one of the main issues of such models. 

The model derived in [5] consists of a Cahn-Hilliard system with transport and reaction terms 
which governs various types of cell concentrations. The reaction terms depend on the nutrient con¬ 
centration (e.g., oxygen) which obeys to a quasi-static advection-reaction-diffusion equation coupled 
to the Cahn-Hilliard equations. The cell velocities satisfy a generalized Darcy’s law where, besides 
the pressure gradient, appears also the so-called Korteweg force due to the cell concentration. 

Numerical simulations of diffuse-interface models for tumor growth have been carried out in 
several papers (see, for instance, O Chap. 8] and references therein). However, a rigorous mathemat¬ 
ical analysis of the resulting PDFs is still in its beginning. To the best of our knowledge, the first 
related papers are concerned with a simplified model, the so-called Cahn-Hilliard-Hele-Shaw system 
(see [13], cf. also [T31 [H]) in which the nutrient n, the source of tumor St and the fraction So oi 
the dead cells are neglected. Moreover, very recent contributions (see |ll[ini[Sll3]) are devoted to the 
analysis of a newly proposed simpler model in m (see also m)- In this model, velocities are set to 
zero and the state variables are reduced to the tumor cell fraction and the nutrient-rich extracellular 
water fraction. 

In what follows we briefly introduce the model proposed in [3], where a complete description 
as well as numerical simulations are provided. Our multi-species tumor model includes the mechanical 
interaction between different species. The following notation will be used: 

• = 1,2,3: the volume fractions of the cells: (pi = P: proliferating cell fraetion; (p2 = (po- 

dead cell fraetion; = pH ^ host cell fraction; 

• H: the cell-to-cell pressure; 

• u:=Ui,i = 1,2,3: the tissue veloeity field. We assume that the cells are tightly packed and they 
march together, 

• n: the nutrient concentration; 

• $ = pT) + P'. the volume fraction of the tumor cells which is split into the sum of the dead tumor 
cells and of the proliferating cells; 

• Jp. the fluxes that account for mechanical interactions among the species; 

• Si,i = 1,2,3: account for inter-component mass exchange as well as gains due to proliferation 
of cells and loss due to cell death. 

The variables above are naturally constrained by the relation pn -|- $ = 1. 

The volume fractions obey the mass conservation (advection-reaction-diffusion) equations: 


dtpi diVa,(U(/)i) = -diVa;Ji -I- 


( 1 . 1 ) 
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We have assumed that the densities of the components are matched. Notice that unlike in [3], for 
simplicity, the variable 4>w standing for the volume fraction of water has been omitted. The total 
energy adhesion, supposed independent of ipH, has the form 

= + dx, 

where is a logarithmic type mixing potential (cf. (12.411 in Subsection l2.1|) . Then, we define the fluxes 
J$ and Jif as follows: 

+ J 2 := -V, = -V, - A$) := -V,/i, 



where we have used in the last equality the fact that (pH = 1 — $ and where /r is the chemical potential 
of the system. For the source of mass in the host tissue we have the following relations: 

• St = Sd + Sp := ^2 + 5i, 

• ^Sh := ^83 = (PhSt = (1 ~ 4?)S'7’. 

Assuming the mobility of the system to be constant, then the tumor volume fraction and the host 
tissue volume fraction (j)H obey the following mass conservation equations (cf. (11.11) '): 

dt^ + div:j(u$) = -diva,J$ + $(S'2 + Si), (1.2) 

dt(l>H + divj,(u(()//) = -div2;J// + $ 53 . (1.3) 

Using now the fact that St = Si + S 2 and recalling that (jin + 4) = 1, we can forget of the equation 
for (j)H and we recover the equation for $ in the form 

+ div,;,(u$) - divj,(VxM) = ^St, M = -F'($) - A$. (1.4) 

As in m, we suppose the net source of tumor cells St to be given by 


St = ST{n, P, 4>) = XMnP — Al(4> — P), 

where Am > 0 is the mitotic rate and Xp > 0 is the lysing rate of dead cells. The volume fraction of 
dead tumor cells (pn would satisfy an equation similar to dni), namely 

dtpD + div2,(u0D) - = ^Sd, 

where the source of dead cells is taken as 

Sd = Sain, P, 4)) = (Aa + XNH{nN — n)) P — Al(4) — P). 

However, we prefer to couple the equation for 4> with the one for P = 4> — pp which then reads 

dtP + diva;(uP) = $(S't - Sd)- 

Here XaP describes the death of cells due to apoptosis (cf. [5J p. 730]) with rate Aa > 0 and the term 
XnHItin — n)P models the death of cells due to necrosis with rate Xn > 0. In [15] H was originally 
taken as the Heaviside function. Here, for mathematical reasons, we smooth it out by taking it as a 
regular and nonnegative function of n. The term un represents the necrotic limit, at which the tumor 
tissue dies due to lack of nutrients. 

The tumor velocity field u (given by the mass-averaged velocity of all the components) is 
assumed to fulfill Darcy’s law: 

u = -VajH-I- 
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where, for simplicity, the motility has been taken constant and equal to 1. Summing up equations 
(HU, we end up with the following constraint for the velocity field: 

diva;U = St- 

Since the time scale for nutrient diffusion is much faster than the rate of cell proliferation, the nutrient 
is assumed to evolve quasi-statically: 


—An + vijnP = Tc{n^ $), 
where the nutrient capillarity term is 

Tc(n, $) = [vi{l - Q{<^)) + V 2 Q{<^)] (ric - n), 

vu represents the nutrient uptake rate by the viable tumor cells, z/i, z /2 denote the nutrient transfer 
rates for preexisting vascularization in the tumor and host domains, and ric is the nutrient level of 
capillaries. The function Q(d>) is assumed to be regular and to satisfy vi{l — (5($)) + y 2 Q{^) > 0 
(cf. (13.21) below). 

Remark 1.1. We chose the boundary conditions proposed in m for /i, If and n. On the other 
hand, under the homogeneous Neumann boundary conditions suggested in m for P, we could not show 
that the system is well-posed. For this reason, we ehose the boundary conditions which are 

natural in connection with the transport eguation for P. In particular, the proliferation function 
at the boundary has to be nonnegative on the set where the velocity u satisfies \i v > 0, with v denoting 
the outer normal unit vector to the boundary of our domain fl. By maximum principle, this implies in 
particular that P > 0 in II, which is an information we need for proving well-posedness of the system. 


In summary, let C be a bounded domain and T > 0 the final time of the process. For simplicity, 
choose Am = vu = I, Xa = Ai, Ajv = A 2 , Xl = A 3 . Then, in fl x (0,T), we have the following system 
of equations: 


+ div3,(u$) - diVx(Va:At) = p = 

(1.5) 

U = -Vxn + /iVa:$, 

( 1 . 6 ) 

II 

(1.7) 

dtP + div 3 ;(uP) = 4 >(S't - Sd), 

( 1 . 8 ) 

— An + nP = Tc{n, $), 

(1.9) 

where 

5T(n,P,$) = nP-A3($-P), 

( 1 . 10 ) 

SDin,P,^) = (Ai+A2P(n7v-n))P-A3($-P), 


Tc{n, 4>) = [ 1 / 1(1 - (5($)) + n 2 Q{^)] {uc - n). 

( 1 . 11 ) 

The functions Q, H and the constants Xi,Vi will be described in Section [2T] System (11.5111.91) will be 
coupled with the following boundary conditions on dO, x ( 0 ,r): 

/i = n = 0 , n = 1 , 

( 1 . 12 ) 

Va,® ■ n = 0, Pu ■ n > 0, 

(1.13) 

and with the initial conditions 

4>(0) = $ 0 , ^*(0) = Po in n. 

(1.14) 

Note that, as P > 0. the second condition in (11.131) should be interoreted as P = 0 whenever vl-v < 0. 
meaning on the part of the inflow part of the boundary. Moreover, in the weak formulation, that 
condition will be incoroorated into eouation (11.81) turninef it into a variational ineonalitv fcf. (12.181) 
below). 

The different nature of the four equations as well as their nonlinear coupling (especially due to 


the Korteweg term in the pressure equation) make the analysis of the problem particularly challenging. 
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Moreover, we may notice that the singular limit studied in the last Section [5] as the interface energy 
coefficient is let tend to zero can be obtained only under more restrictive assumptions on the potential 
which is required to be strictly convex, and under different boundary conditions for u (namely, we 
assume no-flux, rather than Dirichlet, conditions for 11). We refer the reader to Remark l6.II below for 
further comments and for the discussion of related open problems. 

Plan of the paper. The main results and assumptions are stated in Section [5] The subsequent 
Sections [3] and 0] are the core of the paper where we provide the a priori bounds for our solutions and 
we show the weak sequential stability properties. In Section [SJ we construct an approximation scheme 
compatible with the apriori estimates and prove its well-posedness. In the last section, we analyze the 
singular limit problem mentioned above. 


2 Assumptions and main results 


2.1 Singular potential and initial data 

We suppose that the potential T supports the natural bounds 

0 < $(t, x) < I. 

To this end, we take T = C + where B G C^(]R) and 

C : K M- [0, oo] convex, lower-semi continuous, C(<I>) = oo for <I> < 0 or $ > I. 
Moreover, we ask that 


CgCVO,!), lim C'($)= lim C'($) = oo. 
$-> 0 + 

A typical example of such C is the logarithmic potential 

<i>log($) -I- (1 — <&) log(l — !>) for $ G [0,1], 

c($) = ■ 


oo otherwise. 


( 2 . 1 ) 


( 2 . 2 ) 

(2.3) 


(2.4) 


Remark 2.1. Condition (12.31) has mainly a technical character and is assumed just for the purpose 
of constructing a not too complicated approximation scheme (cf. also Remark 15.111 . At the price of 
some additional technical work it could be avoided. One may, for instance, consider the case where 
C(<I>) = 7[op]($) (the indicator function of [0,1]^, which does not satisfy (12.31) . 

Regarding the functions Q and H and the constants Xi, Vi appearing in the definitions of St and So-, 
we assume Q,H £ C^(]R) together with 

> 0 for t = 1,2,3, H>0. (2.5) 

[i/i(l-g($)) + j/2Q($)] >0, 0<ne<l. (2.6) 

Finally, we suppose 12 be a bounded domain with smooth boundary in and impose the following 
conditions on the initial data: 

<i>o£H\n), 0<$o<l, C($o) e (2.7) 

Po G L^{n), 0 < Po < 1 a.e. in SI. (2.8) 


2.2 Main result 

Before stating the main result, let us introduce a suitable weak formulation of the problem. We say 
that ($, u, P, n) is a weak solution to problem (II. 5111.14)) in (0, T) x 12 if 
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(i) these functions belong to the regularity class: 

$ G ^“([O, T];H^{il)) n L2 (o, T; (2.9) 

C($) G L°“(0,T;L^(r2)), hence, in particular, 0 < $ < 1 a.a. in (0,T) x fl; (2-10) 

u G L^((0,T) X n;M^), div u G L°°((0,r) x f2); (2.11) 

( 2 . 12 ) 


P G L“((0,T) X 11), 0 < P < 1 a.a. in (0,T) x 11; 
n G L^{0,T;W^’^in)), 0 < n < 1 a.a. in (0,T) x H; 

(ii) the following integral identities hold: 

f ( + ‘hu • VxV? + da; dl = — f $o‘/^(0, •) da: 

J 0 ^ J vl 

for any ip G C“([0,T) x 12), where 

p = -A$ + P'($), u = -V^n + /rVx$, 
div 3 ;U = St a.a. in (0, T) x H; Va,$ • v\gQ, = 0; 

/ f [Pdtp + Pu ■ VxP + $(S't - Sd)p] dx dt> - ( Po(p{d, ■) da; 
Jo Jn Jn 

for any p G C))”([0,r) x H), p\QTt > 0; 

— An + nP = Tc{n, $) a.a. in (0, T) x 12; n\gn = 1. 


(2.13) 

(2.14) 


(2.15) 

(2.16) 

(2.17) 

(2.18) 

(2.19) 


Now, we are able to state the main result of the present paper: 

Theorem 2.1. Let T > 0 be given. Under the assumptions stated in Subsection \2. 1[ the variational 
formulation !i2.15\f2A^) of the initial-boundary value problem admits at least one solution 

in the regularity class 

Remark 2.2. It is worth observing once more that the second boundary condition (11.131) is now 
incorporated into the variational ineguality (12.181) . 


3 A priori bounds 

In this section we establish several formal a priori estimates for our solution. The procedure turns out 
to be rigorous when (smoother) solutions of the approximated problem (l5.391l5T43l) are considered. In 
particular, this happens for the regularized solution constructed in Section below. In this section 
we refer to system (ll.5Hl.14l) and not to the weak formulation (I2.15H2.1^ because actually the a- 
priori estimates should be performed on the regularized problem (I5.391I5T^51) whose solutions are more 
regular than the ones obtained at the limit. 

We start with noticing that, as a direct consequence of our choice of the potential P, the phase 
field function 4) satisfies dZI]). 

3.1 Lower bound for P 

The density function P satisfies the transport equation p. 8 l) . which can be equivalently rewritten in 
the form 

dtP + \i-VxP = -PST + ^{ST-SD) (3.1) 

= P [—-St + 4> (n — (Ai + X2H{nN — «)))] • 

Thus, provided 

P(0, ■) = Po > 0, and P{t,x) > 0 for x G 9f2, u • z/ < 0, 
we can deduce by maximum principle arguments that 

P > 0. 
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3.2 Positivity and upper bound for n 

In order to obtain positivity of n we need 

-nF + Tc{n,(p) = -nP + [vi(l - Q($)) + i^ 2 Q(^’)] {ric - n) 

to be positive (non-negative) whenever n < 0 ; actually, this follows from the hypothesis (cf. (12.61) in 
Subsection \TT\i 

- Q($)) + >0, 0 < n, < 1. (3.2) 

This assumption also implies that n < 1, so we may conclude that 

0 < n(t, x) < 1. (3-3) 

3.3 Upper bound for P 

Since 0 < $ < 1 and 0 < n < 1, by the assumptions provided in Subsection 12 . 1 1 we have 

—$ (Ai -I- X 2 H{nN — n)) < 0. 

Hence evaluating the expression on the right-hand side of (EH) for P = 1 yields 

P [—S't -I- $ (n — (Ai -I- \ 2 p[{nN — n)))] < A 3 ($ — 1) + n{^ — 1). 

Consequently, provided 

0 < P(0, ■) = Po < 1, and 0 < P{t, x) < 1 for a: G Oil, u ■ ;/ < 0, 

it follows that 

0<P(t,x)<l. (3.4) 


3.4 Estimates for the Cahn-Hilliard equation 

The standard estimates are obtained via multiplication of m by /r: 

4-/ +-^(4’) dx-t / dx = - / u • dx, 

Jn J 

where, by virtue of (EH, 

— / u • dx = — |up dx+ HdivjjU dx = — |up dx+ HS't dx. 

Jo, Jq Jq, Jq Jq, 

Consequently, EH reads 

4 / J'($) dx+ j [|V,/r|2 + |u|2] dx= [ n^rdx, 

ut 7n J Jn Jo. 


where 


HS't dx 


< ||S'T||L“(n)l|n||Li(n)- 


Seeing that H solves the Dirichlet problem 

-AH = St - divx{fJ-Px^), n|an = 0, 

we deduce that 

l|n(t, Olli/qn) < II'S't(C Olli^cn) + llA^Va;$||i2(Q.R3), 

where, by means of Gagliardo-Nirenberg interpolation inequality, 

||Wa;4>||i2(Q.R3) < ||A^(t, •)llL*(n) ||Vx4>||i4(Q.R3) 


(3.5) 


(3.6) 
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<c||A.(i,.)||L4(n)||$(t,-)lli'2(n)l|A$(t,-)lli((a) 


< c||M(t,.)IU4(a)||<i>(i,-)liyi(o) {M%]n) + m\\%]n)) > 


where the last inequality has been obtained testing the second (nm by $ and using the properties 
of J- (in particular, the monotonicity of C'). 

Thus, going back to (13.611 and applying a standard version of Gronwall’s lemma, we deduce 
the bounds 

sup ||$||Hi(n) < c, (3.7) 

te(o,T) 



ll^2:Mlli2(f2.K3) + |up 


dt < c. 


(3.8) 


3.4.1 More estimates on d) 

Knowing that 

- Ad> + C'(d>) = 5 = M - 'B'(d>) G L\0, T; H\n)), (3.9) 

we may multiply this relation by —Ad> and use once more the monotonicity of C' to deduce 

/ ll‘5’ll^2,2(n) dt < c. 

Jo 

Next, take an increasing function h and multiply (13.91) by h{C'{^)) to obtain 


[/i'(C'(d>))C"(d>)|V,,d>p + h(C'($))C'(d>)] dx= gh{C'i^)) dx. 


(3.10) 


Choosing h{-) = (•)® and using that g G L^(0, T; L®(n)), we then easily deduce 

C'($) is bounded in L^(0, T; L®(r2)), 
whence, comparing terms in (021), we also infer 


/ II‘J’II^ 2 , 6 (o) dt < c. 
Jo 


(3.11) 


3.4.2 Estimates on u 

Note that we already know 


diva;U = St bounded in L°°{{0,T) x Q) 


and 

u bounded in L^({0,T) x fl;R^). 

Next, we compute 


curl,,u = A G L^{0, T; L^(fl)) n L\0, T; L^{n)). 

Hence, we may take a test function (p G C'°°(K^) with support contained in fl and apply m p- 51] 
to the function t/ju. In view of the fact that div 2 ;((^u) and curl(y)u) are bounded in L^(0, T; L^(]R^)), 
we then obtain that pu is bounded in L^(0, T; i7^(K^)). Consequently, u satisfies 

rT 


0 


(3.12) 




4 Weak sequential stability 


Suppose that 


{$ 5 , U 5 , Ps, ns}s>o 

is a family of solutions complying with the a priori bounds obtained in the last section. Our goal is 
to show the precompactness of this family of solutions, that is to prove that 


{ $ 5 ^-$weakly-(*)inL°°(( 0 ,T)xO), ' 
U 5 —>■ u weakly in L^((0,T) x 0;R^), 
Ps^ P weakly-(*) in L°°((0,r) x O), 
rii ^ n weakly-(*) in L°°((0,T) X O), 


where the limits solve the same system of equations. 


4.1 Compactness of the time derivatives 

It follows from (ll.5|) and the a-priori estimates we have on $ that 

dt^s —>■ weakly in L^(0, T; 

whence, in accordance with (13.111) and the uniform bounds obtained before, we get 

Va;<l >5 —>■ Va;<l> in L'*((0,T) x for a certain q > 2, (4-1) 

and 


$5 —$ a.a. in (0,T) x fl. 


(4.2) 


Consequently, we can pass to the limit in (II.5L using the fact that diVxU ,5 = St,s and the standard 
monotone operator theory to handle the limit in fj,s- 

Let us now test dEHl) hy (j) G Then, integrating by parts and using (13.811 . we easily 

arrive at 

rT 


\Pt 


1 ^- 


dt < c. 


2(n) ^ 


(4.3) 


Coupling this with (El), we infer 

Ps ^ P strongly in L^( 0 ,T : lT“'^’^(n)) for every e e (0,1). 
Let now (p € with support in il. Then, from (14.41) and (13.1211 . we obtain 

pT p pT p 


/ / Psusp dxdt^ / Pup dx dt, 


Jo Jn Jo Jn 

whence we can identify the limit of the product 

U 5 P 5 -G uP weakly in L^((0,T) x f2;]R^). 
Next, testing dm) by ns and using dm and dm, it is easy to infer 

ns ^ n weakly in L^(0, T; WQ^’^(n)), 


(4.4) 


(4.5) 


(4.6) 


whence, using dm again, 


Psns —>• Pn weakly-(*) in L°“((0,T) x 11), (4.7) 

Psb{ns) —>■ Pb{n) weakly-(*) in L“((0,T) x H), (4.8) 

for any function b, where 5(n) denotes a weak limit of {5(n5)}5>o- 
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4.2 Strong convergence of the nutrients 

We finish the proof of compactness by showing strong (a.a.) pointwise convergence of the nutrients 
{?^ 5 } 5 >o- We have 

- Ans + Pns + [viil - (5($)) + i' 2 Qi^)] ng = (4.9) 

= (-P - P5)ns + [Ml - Qi^s)) + V2Q{‘^s)\nc 

+ [z/i(l - (3($)) + V 2 QM] - - < 5 (^ 5 )) + V 2 Q{^s)] 'jns, 

and, for the limit system, 

-An + Pn+ [^ 1(1 - (5($)) + V 2 QM] n = [vi{^ “ 0(‘^’)) + ^ 2 Q{^)\ ric- 

Thus, testing respectively by ng and n, integrating by parts, and making use of the relations (I4.1H4.8I) 
(in particular, (14.81) is exploited with the choice b{ng) = Ug in order to manage the first term on the 
right hand side of (14.91) '). we may show that 

[ f \Vo,ng\'^ + Png + WM-QW) + ^2QM]nj dx 

Jo Jn 

[ [ |Va;np + + [z/i(l - (2($)) + z/ 2Q(4>)] dx dt, 

Jo Jn 

which yields the desired conclusion 

^xng —)• Vxn, ng ^ n in L^((0,T) x fl). (4-10) 


5 Approximation scheme 

In this section we briefly introduce the approximated scheme needed to obtain rigorously the above 
described a priori estimates. This part is quite standard, hence some details are omitted. 

5.1 Local existence by fixed point argument 

Let 

5 e L\{),T-L\n)), p|Ua(o.T;L2(n)) < R 2 , (5.1) 

where R 2 > Q (this value can be chosen arbitrarily; for instance we can take i ?2 = !)■ 

Replace St + ^ 3 ^ with S and solve dEMTZl) locally in time by a fixed point argument. The 
following can be proven: 

Lemma 5.1. Let S be given by 15.11) . Let 6 € (0,1/4) and $95 £ 1 — ^]- Then 


there exists Tq £ (0,T] possibly depending on 5 such that the system 

a. 9(4 - + u • - A/r = 0, b. /r =-A$ + P'(4), (5.2) 

u =-Va;!! + ^Va;4, (5.3) 

- An =-div,(/rV,4)+:S-A34, (5.4) 

coupled with the initial and boundary conditions 

^ = U = 0, V„4 ■v = 0 on 911 X (0, T), (5.5) 

m(0) = 0, 4(0) = 4o5, (5.6) 

has at least one solution (4,/i, 11, u) satisfying the regularity properties 

n G L®(0,T; M{n)) . 
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(5.7) 

(5.8) 

(5.9) 


Proof. Let Tg G (0,T] to be chosen below and let 

¥GL^(0,To;W^’^(n)), JI G L\{0,To) X n), 

with 

ll‘J’llL4(0,To;Vyi.4(o)) + ||7l||L4((o,To)xn) < ^1- 

This in particular implies 

||MVa;<i)||L 2 ((o,To)xn) < Q{Ri)- 

Again, i?i > 0 can be chosen arbitrarily. Here and below Q is a computable positive function, 
monotone increasing in each of its arguments. 

Replace d>, Ji and S in equation (15.41) . which becomes 

- AH = -div,(7lV,l) + ^ - A3$ (5.10) 

and is still endowed with the boundary condition H = 0. Clearly, (15.1011 has one and only one solution 

UGL\Q,To-,Hl{Sl)). (5.11) 


Moreover, 

Set 


l|n||L2(O.To;Hi(a)) < Q{RiiR2)- 


u —VxH + G L^((0,Tg) x H;K.^), 


(5.12) 


and replace it in (15.2L Once u is assigned we can easily prove existence of a solution to (15.211 . Note 
that no regularization of J- is required. The regularity class of the solution can be formally determined 
multiplying (15.21) a. by ^ and (15.21) b. by Note that 



u • dec 


A ||u||L2(n;R3) ||Ai||L4(Q) II Va;4>||i4(Q.R3) 


< c||u||^2(Q.R3)||Ai||ffi(n) + c||‘b||ff2(Q) 

< c||u|||2(n;R3)||Ai||ffi(n) + c||/^||i2(n) + c||4>||^i(f2), 

The last inequality follows by multiplying dEll) b. by and using the monotonicity of C (cf. Sub¬ 
sec. HD). Then, we can apply Gronwall’s lemma to obtain 

ll‘J^I|L“(0.ro;ffi(n)) + ||M||L“(0,To;i?i(n)) < Q{Ri,R2tS~^ ,To). 

Next, multiplying ()5.2p a. by fit, the time derivative of ()5.2I) b. by and summing the results, yields 

ll‘J’llffi(O.To;ffi(n)) + ll/^llffi(O.To;ffi(n)) < Q(^l:-^2, <5 ^,Tg). (5.13) 

Note that, to deduce (15.131) in a rigorous way, it would have been necessary to regularize C in order for 
its second derivative to be well defined. However this is a standard argument and the resulting estimate 
would be independent of the regularization since it just relies on the monotonicity of C. Hence, we 
omit giving details. Finally, the same argument used for the complete system yields (cf. (13.111) 1 

ll‘J’llL“(0.To;W2.e(n)) < Q(^lj ^2, ^ \7d). (5-14) 

Next, multiplying dO]) a. by A/i and using (I5.14|) . it is not difficult to arrive at 

llMllL“(o,To;ff2(n)) A 0(Ri,-R 2 ,<5 ^,To). (5.15) 

By Sobolev’s embedding this implies that there exists Cs > 0 such that 

— Cg < fJ-(t,x) < Cs for a.e. (t,x) G (0,T) x H. (5.16) 
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Thanks to assumption (E3D, recalling that ^ 1“^]? applying maximum principle arguments 

in (ESI), we deduce the following separation property: there exists > 0 such that 

— 1 + Ki < x) < 1 — Ks for a.e. {t, x) € (0, T) x (5-17) 

Thanks to the above separation property, we can easily prove the uniqueness of the couple ($, fx). We 
already know that, once S is given, a unique 11 solving (15.101) is determined. Hence we have a unique 
u given by (15.121) . Assuming that, for this u, a couple of pairs ($,/r) solve (15.2L we can test the 
difference of (15.21) a. by the difference of the /i’s and the difference of (15.21) b. by the difference of the 
<l>t’s. Performing standard manipulations and using the separation property (15.171) it is then easy to 
deduce a contractive estimate. Hence, the couple ($, p) is in fact unique. 

The above argument permits us to define, for the fixed S given by (ED, the map 

Ml : Br, ^ L°°{0,To;W^’^{n)) X L^(0,To;H^{n)), Mi : (I,/!) ^ ($,/r), 

where Br^ is the closed ball of radius i?i in the space L'‘(0, T; x L^((0,T) x H). We aim to 

apply Schauder’s fixed point theorem to the above map. To this purpose, we first observe that we can 
take To small enough so that the map takes values into Br^^. Moreover, is compact by Sobolev’s 
embeddings. Finally, the continuity of Ali can be shown by standard methods relying on the a priori 
estimates obtained above. 

Hence, by Schauder’s theorem, there exists a time Tq < T, possibly depending on S, such 
that system (EHEl, coupled with the initial and boundary conditions, has at least one solution 
($, /r, H, u), in the interval (0, Tq). The regularity of this solution is specified by (j5.11|) . (15.131) . (I5.14L 
and (15.151) . To conclude the proof, it remains to improve the regularity of u. Since we now know that 
H solves (|5.4I) . using (15.141) and (15.151) it is easy to check that 

II — diVa;(/rVa;$)||icx,(Q^7'g.i6(f2)) < Q{R2,S ^,To). 

Hence, recalling ED and applying elliptic regularity to (ED, we arrive at 

l|n||L»(0,To;ff2(0)) + ||u||2,8(o,To;ffi(f2;K’^)) — ^,To). (5.18) 

This gives (15.91) and concludes the proof of Lemma 15.11 □ 

Remark 5.1. As already noted in E.emark \2.1\. assumption (12.31) is needed only for the sake of 
obtaining higher regularity of approximating functions. Indeed, it would be enough to assume it to 
hold in the approximation (for a suitable family Cs tending to C as S ^ 0) and not necessarily for C. 

Lemma 5.2. For any S as in (ED, the guadruple ($,p,, H, u) solving i5.‘A\5.4\ ) with the initial and 
boundary conditions 15. 5115.61] is unique. 

Proof A contractive estimate can be obtained simply by multiplying the difference of (15.21) a. by the 
difference of the p,’s, the difference of (15.21) b. by the difference of the 4>(’s, and the difference of the 
(ED by the difference of the H’s. We leave the details to the reader. We note that the separation 
property (15.171) and the additional regularity (15.181) play a role in this argument. □ 

Thanks to the above Lemmas, given S, there exists a unique quadruple ($,/i,n, u) solving (15.2115.41) . 
We now plug this quadruple into (a proper regularization of) system (ll.8m.9() . Namely, we have the 

Lemma 5.3. Let S as in 15. 1[) and let Tq, 4>, p, H and u be given by Lem,m,a, \5. 1\ Let Pqs G IlQ(fl), 
Pqs G [0,1] a.e. . Then there exists one and only one couple {P,n) satisfying the system 

dtP — SAP + div 2 ,(uP) = $(n — Ai — X 2 H{nN — n))P, (5.19) 

- An + nP = Tcin,^), (5.20) 

over (0,To), together with the initial and boundary conditions specified at the beginning (with Pq 
replaced by Pqs) and the additional condition 

6P = 0 on on. (5.21) 


12 
























Moreover, 


(5.22) 


a. P{t,x) >0, b. 0 < n{t,x) < 1 for a.e. {t,x) G (0,To) x 
and the following regularity properties hold 

\\P\\m{ 0 ,To-,L'^in)) + ||-P|U“( 0 .To;L 2 (n)) + \\P\\L'^{0,To-,H^{n)) < QiR 2 ,S~^,To), (5.23) 

ll’^l|Hi(0,To;Hi(f2)) + ll’^lli”(0,To;i/2(a)) < Q{R2, S~^ ,To), (5-24) 

Proof. Let us introduce the truncation operator P(r) = max{0, min{l, r}}. Plugging 7” into the right 
hand side of (15.191) . we obtain the elliptic-parabolic system 

dtP - SAP + div^(uP) = $(r(n) - Ai - A 2 i7(njv - n))P, (5.25) 

- An + nP = Tc{n,<i>), (5.26) 

Existence of solutions to (the initial-boundary value problem) for (j5.251l5T^Sl) is standard. For instance, 
one may prove it by using the Faedo-Galerkin scheme. Hence, we omit the details. Rather, we point 
out which are the main a priori estimates involved, with the purpose of establishing sufficient regularity 
properties of solutions. We will also see, as a byproduct, that the component n turns out to take values 
in the interval [0,1] so that the couple (n, P) will in fact solve (I5.191lffi20ll (without truncation). 

To carry out this program, we start with multiplying (I5.25P by P to get 

< f Pu-VxPdx+ ( ^{T{n) — Xi — \ 2 H{nfq — nfjP'^ dx. 

Jq Jq 

where we used condition (15.211) . In view of the smoothness of $ and the presence of the truncation 
operator, the only term that needs to be estimated is the first one on the right hand side. By Poincare’s 
and Young’s inequalities, we have 


Pu • VxP dx < ||P||L4(Q)||u||i4(Q.R3)||Va;P||L2(f2) 


<ll^lll/4%)||P|li(tn)l|u||L4(0;«3) 

< + c5||P||i2(f2)||u|||4(Q.R3). 


Hence, by (15.181) and Grbnwall’s Lemma, we arrive at 

l|-f’llL“(0.To;L2(n)) + \\P\\L^{0,To-,H^(n)) ^ Q{R2,S \To). 
Combining this relation with (15.181) we infer 

||diVa;(uP)||^8/5(o,To;L3/2(n)) <Q{R 2 ,S ^,To), 

whence, applying parabolic regularity theory to (I5.25|) . 

ll-Pi||L8/3(0,To;L3/2(o)) -f ||P||L8/5(o,To;W2>3/2(n)) < Q{R2,S ^,To). 

In turn, by interpolation, this gives 

for all e G (0,1), whence, by Sobolev’s embeddings. 


(5.27) 


(5.28) 


(5.29) 


(5.30) 


(5.31) 


l|4^llL8/3(0,To;L6-qn)) + ll^a;7^llL8/3(0,To;L3-qn)) <Q{R2,S ^,To,e ^). 


(5.32) 
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Consequently, by (I5.18I1 . 


\To,e i), 


whence, going back to (I5.25L 

ll^t" 


L2(0,To;Lt“'(n)) )• 


(5.33) 


(5.34) 


We proceed by bootstrapping. Actually, some more iterations (whose details are omitted for brevity) 
permit us to obtain (15.231) . Once sufficient regularity is achieved, the same maximum principle argu¬ 
ment used for the coupled system gives ()5.22ll a. 

We now pass to equation (15.261) . By elliptic regularity (i.e., multiplying by n — 1 € i?g(f2)), 

we infer 

ll’^l|L2(o.To;i/i(f2)) ^ 0(^2, To). (5.35) 

Next, using, as for the complete system, the sign condition on the right hand side, we get the second 
()5.22p . This entails in particular that P solves (I5.19|) . i.e., no truncation in fact occurs. 

By (I5.28|) . (15.221) . and elliptic regularity, it follows that 

ll’^llL“(0,To;if2(o)) < <5(^2, <5 ^,ro)- (5.36) 

Now, we differentiate (15.201) in time. Recalling (11.111) we have 


- Ant + ntP + nPt = - [vi{l - Q(<I))) -|- V2Q{^)] nt (5.37) 

+ [-i^iQ'($) -b j22 Q'(4')] $t(nc - n). 

Test the above relation by nt and use the regularity given by (15.131) and ()5.23l) together with the 
second ()5.22|) and the positivity of the given term [vi{l — Q(4))) -b i/ 2 Q($)], to obtain 

ll’^t||L2(Q < (5(7?2, <5 ^,To). (5.38) 

This, combined with (|5.36L vields (I5.24L Finallv. we have to Drove unioueness of the solution (n,P). 
This in fact follows from a standard argument. Indeed, it is sufficient to test the difference of ()5.19l) 
by the difference of the P’s and the difference of (15.201) by the difference of the n. Then, the transport 
term in (I5.19|) is treated in a way similar to (I5.27|l . whereas the right hand side of (I5.20|) is easily 
controlled in view of the high regularity of $ and of the sign condition. This concludes the proof of 
the lemma. □ 


We can now finalize our fixed point argument for the complete system. 

Theorem 5.1. Let 5 G (0,1/4), $o 5 G $05 G [5,1-5], Pos G Pqs G [0,1] a.e.. 


Then there exists Ti G (0,r] possibly depending on 5 sueh that the system 

- 5A^t-b u • V^4> - = 0, /r =-A4>-bP'(4>), (5.39) 

U = -Va;!!-b (5.40) 

— diva;U = St = nP — A 3 ($ — P), (5-41) 

dtP — SAP + div 2 ;(uP) = <i)(n — Ai — X 2 H{nN — n))P, (5.42) 

- An + nP = Ta{n,<^>), (5.43) 


coupled with the initial and boundary conditions lil.PA^TTT^ (with <i>o 5 and Pqs replacing $0 and Pq) 
and i5.21\) . has at least one solution (<i>,/r, u, P, n) defined over the time interval (0,Pi) and satisfying 
the regularity properties {5. 7|45.<yi] . h5.18\) . h5.23l\^S7S^ . 
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Proof. We let S be as in (EB, where the choice of i ?2 > 0 is in fact arbitrary. Then, applying first 
Lemmas 15.11 15.21 and then Lemma 15.31 we obtain a unique quintuple ($,/i, u, P, n). Thus, we can 
consider the map 

M 2 -.S ^ S := nP + X 3 P. (5.44) 

In view of (|5.23|) . (15.241) and Sobolev’s embeddings, it is easy to check that 

l|5||L~(0,T„;L=(n)) < Q{R 2 , S-\To), (5.45) 

In particular, we can choose Ti £ (0,To] such that S lies in the closed ball Bp>.^ of L®(0, Ti; L^(r2)). 
Moreover, continuity and compactness of the map M 2 in the topology of L®(0, Ti; L^(n)) are an easy 
consequence of the regularity properties (15.231) . (15.241) . the Lions-Aubin theorem, and the a priori 
estimates in Lemmas 15.1115.31 Hence we can apply once more Schauder’s theorem to M 2 , which gives 
that 

S = S = nP +X 3 P in {0,T) xn. (5.46) 

Hence, (I5.3H5.4I1 reduce to (I1.6H1.7L where St is given by (11.101) . This concludes the proof of the 
theorem. □ 


In order to complete the proof of Theorem 12.11 we need now to pass to the limit in the regularized 
system as d \ 0, assuming of course that $05 —>■ $0 and Pqs Po in suitable ways. We just briefly 
comment on the most delicate part of this step, which consists in the passage to the limit in (15.421) 
in order to recover H2.18|l . The other parts are indeed standard since it can be immediately seen that 
the a-priori estimates performed in Section |2] are still valid on the regularized system and they turn 
out to be also independent of S. 

Taking a test function (p as in (12.181) we multiply (15.421) by (p to obtain 



[Pdtp + Pu ■ VxP -I- $ (n - Ai 


X 2 pl{nN — n)) p] da; dt 


(5.47) 


= d 


[ f Va;P • VxP da; dt — / Pop{0, ■) dx + f j Pu • np 
Jo Jq Jn Jo Jan 


dS, 


where, as 


-6 / VxP ■ np dSx, 

Jo Jan 

P\an = d, P>0in (0,r)xO, 


tan 


Pu ■ np dSa; = 0, 




VxP ■ np dSa; > 0. 


Letting 5 \ 0 in (15.471) we get (12.18|1 . 

Finally, let us notice that, by standard arguments, it is possible to show that the a priori 
estimates provide an extension of the local approximate solution up to the original final time T. 
Hence, in particular we have a global solution in the limit. This concludes the proof of Theorem 12.11 


6 Singular limit 

In this section, we consider the problem obtained from (11.5H1.911 by taking St = Sd = d. Hence we 
just consider the system of equations for $ and u, decoupled from the rest, of the form 

-I- div3,(u4>) - diva;(Va:/a) = 0, /a = -e^A$ -|- P'($), (6.1) 

n=-X/xB +pVx^, ( 6 . 2 ) 

diva;U = 0, (6.3) 

with the boundary conditions 

u • n\dn = 0, X/x^ ■ n|an = 0, p\on = 0. (6.4) 
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Notice that, in particular, we are considering here a no-flux condition for If in place of the Dirichlet 
condition in (11.1211 . 

Similarly to Section 13.41 we derive the energy balance 


d 

dt 


|V,$p+.F($) 


dx+ + |up dx = 0. 

JQ 


Next, 


f [e 2 |A$| 2 -f.F"($)|V,$| 2 ] dx= [ V,/x-V,$dx. 
Jn Jn 

Then, assuming strict convexity of T, namely 

T" >X> 0, 


the following estimates can be deduced 

rT 


f dt < c, f ||Va;<l>|||2(Q.R3) dt < c. 

Jo Jo 


(6.5) 


( 6 . 6 ) 

(6.7) 


6.1 Compactness of the velocity 

In view of (16.5p we may assume there is a subsequence such that 

Ug u weakly in L^((0,T) x fl;]R^). 


Obviously, 

We can now write 
whence, seeing that 
we conclude that 


div:j,u = 0, u • n|an = 0. 

W = -V,, (n^ - ^■(I'e)) - 

e^A$eV,r$e ^ 0 in L^((0,T) x ft), 

curhu = 0, 


which, combined with (16.8|) . yields 

u = 0. 


Therefore, taking e ^ 0, system (EH)-(ESI) converges to 

dt^-Aii = 0, /r = .F'($), 


( 6 . 8 ) 


(6.9) 


which satisfies the energy law 

4 / .F($) dx+ [ dx = 0. (6.10) 

dt Jn Jn 

Summarizing, we have proved the 

Theorem 6.1. Let the assumptions given in Subsec. [All hold, let satisfy (16.61) . and let ($£,^E,Ue) 
denote a family of weak solutions to the system complemented with the boundary eonditions 

EH and the Cauchy conditions. Then, os e —>■ 0, the functions (<i>e,/i^, u^) suitably tend to a 
triple ($,p,, 0 ) satisfying (16.9p together with the energy equality 1 ) 6 .lOp and the initial and boundary 
conditions. 


Remark 6.1. It would be interesting to investigate whether similar estimates could be derived for the 
singular flux 

u = -Vj;!! -b 

e 

However, the above argument does not seem to be easily adaptable to cover such a situation. For 
instance, we cannot prove uniform integrahility of the product 

in that case. 
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